LU-DECOMPOSITION OF A NONCOMMUTATIVE LINEAR 
SYSTEM AND JACOBI POLYNOMIALS 



OSCAR BREGA AND LEANDRO CAGLIERO 

Abstract. In this paper we obtain the LU-decomposition of a noncom- 
mutative linear system of equations that, in the rank one case, character- 
izes the image of the Lepowsky homomorphism U(g) K — > U(t) M ® U(a). 
This LU-decomposition can be transformed into very simple matrix iden- 
tities, where the entries of the matrices involved belong to a special class 
of Jacobi polynomials. In particular, each entry of the L part of the 
original system is expressed in terms of a single ultraspherical Jacobi 
polynomial. In turns, these matrix identities yield a biorthogonality 
relation between the ultraspherical Jacobi polynomials. 



1. Introduction 

1.1. The noncommutative linear system. Let k be a field of character- 
istic zero, let A be an associative, not necessarily commutative, k-algebra 
with unit and let M. be a unital ^4-bimodule. 

This paper is devoted to perform the gaussian elimination process on the 
following non commutative homogeneous system of infinitely many linear 
equations and infinite unknowns 

Ep + Epi + Ep 2 ... = PqE - p\E + p 2 E 

E 2 p + 2E 2 Pl + 2 2 E 2 p 2 ... = p E 2 - 2 Pl E 2 + 2 2 p 2 E 2 . . . 



(1.1) 



E 3 p + 3E 3 Pl + 3 2 E 3 p 2 ... = p E 3 - 3p x E 3 + 3 2 p 2 E 3 



where E is a given element of A and Po,pi,p 2 , ■ ■ ■ belong to M and are the 
unknowns. 

It is clear that (po,pi, ■ ■ ■ ,Pd, 0, . . . ) is a solution of the system (jl.ip if 
and only if the polynomial p = Po + pit + p 2 t 2 + • • • + Pdt d : belonging to 
M\t\ := M. <8> k[f], satisfies the equations 

(1.2) E n p(n) = p(-n)E n for all n £ N. 

This system can not be expressed as a single matrix equation AX = 
with A a matrix with coefficients in A. In fact, non commutative systems of 
(homogeneous) linear equations can seldom be expressed as a single matrix 
equation with coefficients in A. Even when it is possible, only in exceptional 
cases a gaussian elimination process can be performed in a satisfactory way. 
For instance, if dim^ A < oo and dim^ M. < oo, by using bases of A and M., 
one can express every finite system of linear equations as a single matrix 
equation AX = B with coefficients in k. However the size of A might be so 
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big that this approach does not help very much. Some papers dealing with 
this subject are jCj. [US]. [Or] or [GGRWj . 

On the other hand, by using the left and right regular actions L, R : A — > 
Endk(A^) of A on M, it is indeed possible to express every non commutative 
system of linear equations as a single matrix equation but with coefficients 
in Endk(-M). In our case, the system (ll.ip can be expressed as single matrix 



equation MnX = 0, where X 




and 



■ Le~Re Le+Re 



-Re 



Ma 



L 2 E -R% 2L 2 E +2R 2 E 4L|-4R| 



L%,-R%, 3L 3 E +3R% 9L%-9R% 



Le+Re 



HL 2 E -\-HR E 



27L|+27R| 



V : : : \ ) 

In order to have an insight into the LU-decomposition of this system assume 
for a moment that the unknown polynomial has degree 2. In this case it is 
easy to see that 



/ Le~-Re Le+Re 



-Re \ 



L E~R-E 



2L 2 E +2R 2 E 4L|-4i?| 



\L%-R% 3L%+3R 2 E 
/ 1 



9L%-9R%J 




Le—Re Le+Re Le—Re 
{L e -Re? 3L 2 e -3R 2 e 











2(L E -R E r 



= Le+Re 1 

\ L 2 e +ReLe+R 2 e 2 Le+2 Re 1, 

is the LU-decomposition of the system The goal of this paper is to 

obtain a transparent expression for the LU-decomposition of the following 
more general system 

(1.3) E n p(H + n) = p{H - n)E n for all n G N, 

for some H £ A. 



1.2. Relevance of the system (II. 3|) . The interest on the system (jl.3p 
comes from the fact that their solution set is closely related to invariant 
spaces under group actions. This is clear in the particular case in which 
H = and E is the identity of A, since the solutions of (|1.3p are the even 
polynomials. It is not difficult to see that in other more general situations 
there exists a group Ge associated to E, acting on M. <g>k k[i], such that 

{p G M ® k k[t] : E n p{n) = p(-n)E n for all n G N} = (M ® k k[t]) Gfi . 

A second, and more relevant, example of the relationship between the sys- 
tem (|1.3|) and invariants of groups appears in a problem from representation 
theory of Lie groups. Let G a be a connected, noncompact, real semisimple 
Lie group with finite center, let g be the complexification of the Lie alge- 
bra of G and let K Q denote a maximal compact subgroup of G Q . By the 
fundamental work of Harish-Chandra it is known that many deep questions 
concerning the infinite dimensional representation theory of G Q reduce to 
questions about the structure of the algebra U(g) K °, the i^ -invariants in 
the universal enveloping algebra of g. In [T], Tirao proved that the elements 
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of the image of U(q) Ko by the Lepowsky homomorphism (see [L]) satisfy a 
system of linear equations completely analogous to system (jl.3p . Moreover, 
for rank one classical Lie groups it is proved in |BCT| that this image coin- 
cides with the solution set of the system introduced by Tirao. In order to 
prove this result, the LU-decomposition obtained in this paper constitute a 
very useful tool. These connections with groups invariants are not treated 
in this paper. 

1.3. Main results. In $2] we find a periodic sequence of elementary row 
operations that, when applied to the system (|1.3p . transforms it into a tri- 
angular system that involves the adjoint action of E and some combinatoric 
sums of discrete derivatives of the polynomial p. The main result of this 
section is Theorem 12.21 As a corollary we obtain that if p G M. [t] satisfies 

E n p(H + n) = p(H — n)E n for < n < deg(p), 



then p satisfies (II. 3p . 

Note that M$ is an infinite matrix with coefficients in Endk(-M) and it is 
closely related to the matrix 



/ x— 1 x+1 x—1 x+1 \ 

c 2 -l 2x 2 +2 4z 2 -4 8x 2 +8 

c 3 -l 3x 3 +3 9x 3 -9 27x 3 +27 

V / 



€ k[x\. 



Taking into account the relationship between Mq and Mq, the results ob- 
tained in Theorem 12.21 can be translated in terms of polynomials in one 
variable with coefficients in k. The combinatorial aspects of Theorem 12.21 
are transformed into very simple matrix identities, where the entries of the 
matrices involved are a special class of Jacobi polynomials. In turns, these 
matrix identities yield biorthogonality relations between the ultraspherical 
polynomials. It has recently appeared in the literature some other matrix 
identities (in particular LU decompositions) involving Jacobi polynomials 
that translate sophisticated polynomials identities into very simple matrix 
identities (see for instance |KO] ) . 

These results are described in what follows. 

Let denote the Jacobi polynomial of degree n associated to the 

numbers a and (3. Usually a and (5 are complex numbers but, since we 
are working in the context of an arbitrary field of characteristic zero, they 
are assumed to be rational. Recall that the polynomials corresponding to 
parameters a = (3 are known as the ultraspherical or Gegenbauer's polyno- 
mials. Let 



P ^(x) = (x- i)»p«e (f±i) 



It is clear that pn (x) is again a polynomial and it can be expressed in terms 
of the hypergeometric function of Gauss as follows (see (4.22.1) in [Sz]), 

p a, P{x) = (n+a+f3+iy..(2n+a+f3) ^ ( _^ _ n _ q; _ 2n _ q _ ^ % _ ^ 

Our results in terms of Mq are the following: 
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(1) The LU-decomposition of Mq is 



M n 



L U 



( 1 



21 



-1 
1 



1 

-15 
25 
-10 
1 



where 



(Lo)ij 
(Uoh 



i > j; 



and the explicit numerical matrices are respectively formed by the fac- 
torial and (up to the minus signs) Stirling numbers of the second kind. 
(2) The inverse of the matrix Lq is 

i'V), ( 1 )' J j !'!■{,{■>■ U i>j- 

This yields the following "discrete orthogonality" relationship that in- 
volves once many of the ultraspherical Jacobi polynomials with integer 
parameters 



i pi, 

1M) 
2*1 



3*2 



2*0 



V 



1,1 2 p2,2 3 „3,3 
1 3*0 
,1,1 2 p2,2 3 3,3 j 4,4 



4 3 4 2 



5 pa 
4*1 



*i 



t*v •/ 



2,-2 



-2,-2 



■3, — 3 rj — 3, — 3 rj — 3, — 3 
*1 *0 



p 3 - 4 '- 4 p- 4 '- 4 p- 4 '- 4 



7 



V 



A 
7 



(3) We mentioned above that we performed to the system (II. 3p a gaussian 
elimination process in which we periodically apply the same sequence of 
elementary row operations. This periodicity is translated in the follow- 
ing terms. Given an infinite matrix 



an ai2 ai3 

(121 0-22 123 
»31 a 32 0,33 



let s(T) be the shifted matrix of T, that is s(T) is the diagonal blocked 
matrix ( 1 T ) formed by the identity matrix of size lxl and T. We say 
that an infinite matrix A admits a periodic gaussian elimination process 
if there exist an lower triangular matrix with l's in the diagonal To such 
that the sequence of matrices 

A, T A, s(T )T A, s 2 (To)s(To)T A, ••• 

converged to an upper triangular matrix. In this case To is called a 
fundamental sequence of the Gaussian elimination process of A. It is 
not difficult to see that an infinite matrix A admits a periodic gaussian 
elimination process if and only if A admits an LU-decomposition. In 



The concept of convergence is used in the discrete sense, that is a sequence of matrices 
Bfc converges to B if for every i, j there exists ko such that (-Bfc)ij = Bij for all k > kg. 
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fact, if A = LU then the fundamental period is To = s(L)L 1 and, 
conversely, if To is the fundamental period of A then 



V 



s 2 (T ) S (T )T and L = T^ l s(T^ l )s 2 (T Q l ) ■ ■ ■ . 



It is worth mentioning at this point that the Vandermonde matrix 
Vij = V~ x admits a periodic gaussian elimination process and its funda- 



mental sequence is 




and that the recursive definition of 





the Pascal numbers is equivalent to the identity 
= • • • s 2 (T )s(T )T , with T 



It turns out that Mq admits a periodic gaussian elimination process 
and we show that its fundamental sequence is 

/ -x-l 1 • \ 

To = x —x—1 1 ■ I . 
V x -x-l 1 • J 

It is not difficult to see that the inverse of To is the matrix formed by 
the cyclotomic polynomials 



i 

l+x 
1+x+x 2 



1 

l+x 



1 



l+x+x 2 +x 3 1+x+x 2 l+x 1 



and therefore the ultraspherical polynomials can expressed as 



/ i 1,1 
/ jPo 



1 1,1 2 2,2 

2?V 2^0 

1 1,1 2 2,2 3 3,3 

\v 2 fp 

1^1,1 2 ^ 2 ' 2 3 ^ 3 ' 3 l^ 4 ' 4 



I'':! 



4P2 4P1 



Pa 



f 1 



1 

-x-l 1 



\( 1 1 

—x—1 



-x-l 1 
X —x — 1 1 



-x-l 1 

x —x—1 1 ■ 



and 



(Pa 



1,-1 



\ 



-2,-2 -2,-2 
Pi Pa 

-3,-3 -3,-3 -3,-3 
P2 Pi Po 

-4,-4 -4,-4 -4,-4 -4,-4 
P3 P2 Pi Po J 



1 

l+x 1 
1+x+x 2 l+x 1 
l+x+x 2 +x 3 1+x+x 2 l+x 1 



r 
\ . 



1 

l+x 1 
1+x+x 2 l+x 1 



7 



V 



1 

l+x 1 



7 



We notice that in these two identities the polynomials involved are pn' a 
instead of P"' Q 
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2. The gaussian elimination for the system (jl.3p 

Let k be a field of characteristic zero, let A be an associative, not neces- 
sarily commutative, k-algebra with unit and let Ai be a unital ^4-bimodule. 
Given an element r G A, let L r ,R r G Endk(-M) respectively denote the left 
and right actions by r, and let ad r = L r — R r be the adjoint action of r in 
M. 

Recall that M[t] = A4 ®k k[i], and that given an element r € A one has 
the evaluation map M[t] — > M defined using the right action of A on A4 
by a ® p(t) i — ^ ap(r). If r G *4 and p € A4[t] then the evaluation of p in r is 
given by p(r) = po + Pir + • • • + PrJ n ■ 

Given two arbitrary elements E and H in A we are interested in the set 
of polynomials p G M [t] that satisfy the following system of linear equations 
(|1.3p . that is polynomials p G A4[t] satisfying 

E n p(H + n) = p(H - n)E n for all n G N. 

In this section, we assume that E is an eigenvector of ad# in A with eigen- 
value c and we shall perform a gaussian elimination process to obtain, in 
Theorem 12. 2i a triangular linear system in A4 [t] equivalent to the system 

dLU). 

2.1. A first look at the system. We first discuss the linear system (j 1 . 3 j) 
for polynomials of degree zero and one. 

Degree zero: Let p = po G M\t\ be a constant polynomial. Observe that p 
satisfies the first equation of the system if and only if po commutes with E, 
and therefore p satisfies all the other equations. Thus p satisfies the system 
(|1.3p if and only if po G keradg. 

Degree one: A linear polynomial p = po + pit G M[t] satisfies the system 
(|1.3p if and only if the vector (po,pi) is a solution of the following system, 

(2.1) E Po +E Pl (H+l) = p E + Pl (H-l)E 

(2.2) E 2 p + E 2 pi(H+2) = p E 2 + Pl (H-2)E 2 

(2.3) E 3 p + E 3 Pl (H+3) = p E 3 + Pl (H-3)E 3 



It is not difficult to see that equation (|2.3[) can be obtained from the first 
two as follows, 

2E x (O - e 2 x (EB -fix flO x £ - (EH) x fi 2 + (E2J) x e. 



In fact, we shall see in Corollary 12.31 that for all n > 3, the nth equation can 
be obtained from equations (|2.1|) and (|2.2p . Then, in this case, the system 
(jl.3|) is equivalent to the 2x2 system formed by equations (|2.1|) and (|2.2|) . 
This 2x2 system is represented by the matrix 



M 



Le—Re LeRh+i—ReRh-i\ 
yL E —R E L e R h+ 2~R e Rh-2 J 
with coefficients in Endk(TW). Now, since 

Le—Re LeRh+i—ReRh-\ \ I Le—Re LeRh+i—ReRh-i 
L E —R E L E R H+2 — R e Rh-2I \ (Le—Re) 2 
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we obtain that equations (|2.ip and (|2.2p . are equivalent to the triangular 
system U ( b b °) = 0, with 



U 



Le-Re L e Rh+i-ReRh-i\ / Le-Re Le+Re \ f I R h 
{Le-Re? J \ (L E ~R E ) 2 ) \o 1 



In Theorem 12.21 we shall extend this triangularization process to polyno- 
mials of arbitrary degree. In particular, we shall obtain in Corollary 12.31 that 
for polynomials of degree N the whole system f j 1 . 3 j) is equivalent to the first 
N + 1 equations of the system (|1.3|) . 

2.2. Some additional notation. For h G k, h ^ 0, we define the /i-discrete 
derivative % : M [t] — > [£] by 

W) = P(t) "f for P eM[t]. 



It is clear that 



and, if p = p m t m + • • ■ + Po, one has 

d n hP (t) ■ 



0, if n > m 

m\p m , if n = m. 



If (t) n is the Pochhammer polynomial of degree n, that is 

(t) = 1 and (t) n =t(t + l)...(t + n-l) for n > 1, 
it is easy to see that dh(t/h) n = n(t//i) n _i and thus 



(2.4) = (n-r + l) r (t/h) n -r for r < n. 
Given r£iwe extend to M[t] the adjoint action of r in by 

ad r (p + H l-p n i n ) = ad r (p ) + ad r (pi)t H h ad r (p n )i n . 

Observe that ad r commute with the discrete derivative dh for all h € k. In 
addition, if E, H G .4 are elements that satisfy ad# (£7) = ci?, c E k, then it 
is straightforward to prove that 

(2.5) ads (p(i + H)) = adjs(p)(i + ff) - c9 c p(t + 

(2.6) £ n p(t + H) = ("I adjB~ J '(P)(* ~ c i + for all n G N 

(2.7) = (H — nc) m E n for all m, n G N. 
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2.3. The gaussian elimination. For k > and n > k we introduce the 
following equation, 

k n—k 

«£ ■■ E EC" 1 )'® C7*) ad ^°' +i) tfflU P(H + " - cj - i)^ 

i=0 j=0 

= E(~ 1 )* (?) ad £T i a i a c-i p(H-n + 2k- iJ^-C*-*). 
i=0 

We now collect some properties of these equations. 

(1) The equation eq° is, £™ =0 (") ad n f j p(H + n- cj)E^ = p{H - n)E n . 

Hence, in view of the identity (|2.6p . the system {eq° : n € N} is the 
original system (jl.3p . 

(2) On the other hand, the equation eq™ is 

n 

E(-ir(") ad^flp^ lP (H + n-»)^ 

t=0 

n 

= E^" 1 ) 4 (?) ad ^ - n + 2n - 

t=0 

which is trivial for all n £ N. 

(3) The fundamental property of eq^, for k > and n > fc, is that the 
polynomial p appears under a derivative of order grater than or equal 
to k. This implies that for all n > k the equation eq^ involves only the 
coefficients pj with j > k. 

For k > 0, let Eq fe be the system 

Eq fc = {eq?, eq^, . . . , eq^ 1 } U {eq£ : n > k}. 

As we have just observed, Eq° is the original system (II. 3p . The following 
picture may help to keep these definitions in mind. 

k=0 k=l k=2 k=3 k=A k=5 
n=0 o 

n=l o 

n=2 • o 

n=3 • ■ o 

n=4 • ■ o 

n=5 • ■ o 

:::::: > \ equations eq™ 
| f are trivial 

Eq° Eq 3 

Finally, let Eq°° be the system 

Eq°° = {eq™ +1 : n > 0}. 
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According to item (3) above, Eq°° is an upper triangular system, in the 
sense that if p = Ylj>oPj & then eq™ +1 only involves the unknowns pj with 
j > n. The equations of this system are the following, 

n 



eq™ +1 : ^(-l)'® ad^ 1 "' d^U p{H + n+l- i)E* 

n 

+ ^(-1)^") ad^ &tdi_ x p(H + n + l-c- i)E l+1 

n 

= E(" 1 ) i (i) d ? d c~i p(H + n-l- i)E? 



i=0 

n 

Remark 2.1. If deg(p) = N then the only non zero equations of Eq°° are 
eq" +1 for n = 0, . . . , N. In this case, the last two non zero equations are, 

eq# +1 : ^ N e + \pn) = 0, 

eqj^ 1 : adf (pjv-i) + iVadf (p N )(H+^) + N(N+1-Nc) a d£- % (p N )E = 0. 
Indeed, the equation eq^ +1 is 

ad^ +1 Sf p(H+N+l) + adf d? p(H+N+l-c)E = adf p(H+N-l)E, 

however, since p = Nlpjy the equation eq^ +1 becomes adg +1 (pjv) = 0. 

For eq^" 1 we use that d^~ l p = (N-l)\(Np N t - {^Pn + Pn-i), and 
proceed in a similar way. 

The main theorem of this section is the following. 

Theorem 2.2. Assume that [H, E] = cE. For all k > 0, the first N equa- 
tions of the system Eq fc are equivalent to the first N equations of the system 
Eq°. In particular this is true for the upper triangular system Eq°°. 

Proof. The theorem easily follows by induction once we prove that the equa- 
tions eq^ satisfy the following recursion formula 

(2.8) eq£ = eq*" 1 -2 eq^} E + eq^-l E 2 - ad^eq^i) + ad E (eq^)£ 
for n > k > 0. 

The proof of this recursion formula is very technical and will be done 
at the end of this section. We point out that this formula shows that, for 
n > k, the equation eq^ of Eq fc is equal to the equation eq^" 1 of Eq k ~ 1 plus 
a linear combination of the two previous equations of Eq fc_1 . □ 

Corollary 2.3. Assume that [H,E] = cE. If p G M[t] satisfies 
E n p(H + n)= p{H - n)E n for n = 1, . . . , deg(p), 
then E n p(H + n) = p(H — n)E n holds for all n > 0. 

Proof. Let N = deg(p). The assumption on p says that it satisfies the first 
./V equations of the system Eq°, then Theorem 12.21 implies that p satisfies 
the first iV equations of Eq°°. On the other hand, we showed in Remark 12. II 
that p automatically satisfies all the remaining equations of Eq°° . Therefore 
p satisfies all the equations of Eq° . □ 
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In the following corollary we describe the particular cases in which c = 
and c = 1, in this cases the triangularized system became much simpler. 

Corollary 2.4. If ad h(E) = then the equations of the system Eq°° are 



n 

e( ln + i : E [(?) ad ^ ^ P(H +71+1)0 + (?) adr ^ p^+n+l)^ 1 

8=0 

n 

= E (i) ad ^ a ^ v{H+n-l)E* x . 

i=0 



J/ ad// (E 1 ) = E 1 i/ien i/ie equations of the system Eq°° are 



eq™ +1 : ad™ +1 <9? p(ff + n + 1) + &d% c\ n+1 p{H + n)E = 0. 
Proof. If E and H commute then c = 0, hence 



d c -ip(t) = a_u>(i) = + 



Therefore <9*_ 1 p(t) = (— + i) and the formula follows. On the other 
hand, if ad#(i?) = -E then c = 1, and hence <9 c _i = 8q = 0. Therefore all 
the terms of eq™ +1 for i > are zero and thus 



eq™ +1 : ad™ +1 p(H + n + 1) + &d n E d\ p(H + n)E = ad£ p(H + n - 1)E, 



which is equal to the equation stated in the corollary. □ 



2.4. Proof of the Recursion Formula 12.81 We work first on the left 
hand sides (LHS) of the equations involved in the Recursion Formula 12.81 
(see the definition of eq^ in Subsection I2.3[) . In fact, we need to prove that 



LHS( eq^ 1 -2 eq^i E+eq^l E 2 - a d E (e q k n -_\)+ a d E (eqtl)E) = LHS( eq£ ) . 
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Since the sums in the definition of eq^ runs over all possible values of i and 
j for which the binomial numbers are not zero we omit them. 

LHSteojT 1 ) : 

(Y) ("T 1 ) ad T U+l) df'di-i HH + n- cj - i)E^ 1 
LHS(-2eq^ 1 1 S) : 

"("I)* (Y) ("J*) adr dt'OU b(H + n-l-cj- l )E^+ 1 
-(-IT (Y) ("J*) ^7 lHm dt'dU b(H + n-l-cj- l )E^ +1 

LHS( + eqt 1 2 ^ 2 ): 

(n-i-fc) ad 7 2 -^ +i ) fl*" 1 ^ b(H + n-2-cj- i)E^ +2 

LHS(-ad s (eq^ 1 1 )) : 

-(-If (Y) ( n f) <~ St 1 **-! b(H + n-l-cj- i)E^ 

(Y) i n t) dt'dU b(H + n-l-cj-c- i)E^ +1 

+ (-l) i ( fc T 1 ) ( n T fc ) ad^ -1 " '* flf" 1 ^ b(H + n-l-cj- 
LHS( + ad B (eq^)£) : 



+(-ir(Y)r 


-i- 

3 




n n— 1- 

ad E 


-0'+*) gfe- 


l dl 


.! b(H + n-2 


- cj 


- i)E j+i+1 


+(-ir(Y)r 


-1- 

i 




ad s 




- l dl 


_ x b{H + n-2 


- cj 


- c - i)& +i+2 


-(-ir(Y)r 


-l- 

3 




,71-2- 

ad E 




- l dl_ 


_ x b{H + n-2 


- cj 


- i)& +i+2 



We now indicate how to simplify these ten lines. 

Lines 3 and 7 cancel out and so do lines 4 and 10. Lines 6, 8 and 9, after 
a change of variables, become 



"("I)* (Y) ("I?) adr°' +i) dt'dU b(H + n-l-cj- i)E^ 
+ (-1)* (Y) C7Y) ad^ (j+i) dt'dU b(H + n-2-cj + c- i)E^ 
+(-iy (Y) ("7-2") adF dt'dU b(H + n-2-cj + c- i)E^\ 
which is equal to 

"("I) 1 (Y) iT-i) ad E~ U+l) d t ld c-i b(H + n-l-cj- i)E^ 

a d n E ~ U+t) dt'dl^ b(H + n-2-cj + c-i)Ei+ i 

and finally equal to 

(-I)* it 1 ) ( "liO ad^ #YCi &(# + n - 1 - c(j - 1) - (i + 
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Similarly, lines 1, 2 and 5 yield 

(-1)*(Y)C7*) ad^'+^-i KH + n-cj -i)E^ 
"("I) 1 (Y) (T-i) ad ^°' +i) dt'dtt] b(H + n- c(j - 1) - (i + l))E>+\ 
Thus we have, 

("I) 4 (Y) (Jlf) flf" 1 ^ 6(H + n - 1 - c(j - 1) - (i + 

+ (-l)'(Y) (V) ^r (i+i) 6(H + n - cj - i)E^ 

-(-I) 4 (V) {T-D ad F°' +i) ^ _1 S£ti b(H + n- c(j - 1) - (i + 1))^+* 

= "("I) 4 (Y) (-IS adr°' +i) b{H + n- c(j - 1) - (< + 

+ (-irft 1 ) ("7") ad F°' +l) ^-i 6(H + n - cj - i)^* 

= ("I) 4 (Y) ("7") ^-1 b{H + n- cj - i)E^ 

+ ("I)* (Y) ("7") ad E~ ' +i) ^iCi + " " " 

= ("IF © Cj*) a dr °' + ° b(H + n - cj - i)23*+< 

= LHS(eq^). 

We now work on the right hand sides (RHS) of the equations involved in 
the Recursion Formula 12.81 

RHS(eq^ 1 ) : 

(-iy(Y) ad*- 1 -^- 1 ^.! b{H 

RHS(-2eq^ 1 1 S) : 

-(-ir(Y) ad^dt'dU b{H 

RHS( + eqt 2 ^ 2 ) : 

ad^dt'dU b{H 

RHS(-ad £ (eq^: 1 1 )) : 

-(-iy{ k f) ad^T^f" 1 ^! b(H-n + l + 2(k - 1) - i)E n -^ 
-(-ly^T 1 ) adl -1 -*^- 1 ^.! &(# - n + 1 + 2(jfe - 1) - c - 
+(-l) i ( fe 7 1 ) adl -1 -*^- 1 ^! &(# - n + 1 + 2(jfe - 1) - i)f;"-(*-i-0 

RHS( + ad E (eq^)£) : 

ad! - * Sf -1 ^! 6(fT -n + 2 + 2(fc - 1) - i)£™-( fc ^ 
+(-l) i (*7 1 ) ad*- 1 ^^- 1 ^! b{H -n + 2 + 2(k-l)-c-i)E n -( k - 1 ^ 
-(-ly^T 1 ) adl- 1- *^- 1 ^.! &(# - n + 2 + 2(jfe - 1) - fj^-f*" 1 -*) 



- n + 2(fc - 1) - i)^-^- 1 -*) 

- n + 1 + 2(Jfc - 1) - iJ^-C*- 1 -*) 

- n + 1 + 2(fc - 1) - i)^"-**- 1 -*) 

- n + 2 + 2(jfe - 1) - ^E^- 1 -^ 
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Now the sum of lines 1 and 2 is equal to 

-(-l)^; 1 ) ad!" 1 -* d\di_ x b(H -n-l + 2k- i)^-^" 1 -*), 
similarly the sum of lines 5 and 8 is equal to 

(-lJ^V) ad E~ td i d Li b(H-n + 2k- i)£^-C*-0, 
and the sum of lines 6 and 9 is equal to 

(-1)* ({) ad!" 1 -* Sf b{H -n + 2k-c- i)^-^" 1 "*). 
Finally, the sum of the last three lines is equal to 

- ad|~ 1_< b(H -n + 2k-(i + l))^-^" 1 -') 

+ (-^(V) ad l~^i a c-i b{H-n + 2k- i)E n -^ k - 1 
=(-l) < (tD adl - * b{H-n + 2k- i)E n ~^ 

+ (-l)*^ 1 ) ad^af^ b(H-n + 2k- 

=(-l) l '(-) ad|~* SfSj.! b{H-n + 2k- i)E n ~^ 

= RHS(eq^). 
This finishes the proof of the Recursion Formula 12.81 

The following two sections are devoted to express the system (jl.3j) in 
a matrix form and to give its LU-decomposition in a general context. In 
particular we will not assume that [H, E] = cE. 

3. Infinite matrices 
and the matrix representation of the system (|1.3|) 

In this section we recall some preliminaries on infinite matrices, introduce 
the concept of periodic Gaussian elimination process and present the matrix 
representation of the system (11.31) . 

3.1. Infinite matrices. Given an associative k-algebra A, let A4oo{A) de- 
note the set of all infinite matrices 

fa n ai2 ai3 • \ 

. (121 a 22 023 - 

«31 ^32 033 - 

v . . . / 

with aij € .4. The n-minor of a matrix ^4 € .Moo (.A) is the matrix corre- 
sponding to the upper-left corner of A of size n x n. We say that a sequence 
of matrices Bk converges to S if for every i, j there exists ko such that 
{Bkjij = -Bjj for all fc > ko. A matrix A € Moo(A) is said row-finite (resp. 
column-finite) if every row (resp. column) of A contains only a finite number 
of non zero elements. Lower triangular and upper triangular matrices are, 
respectively, examples of row-finite and column-finite matrices. 

It is clear that .Moo (.4) is not a ring since the multiplication of two ma- 
trices does not always exist. Nevertheless, if A, B € M.<x>{A) and either A is 
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row- finite or B is column-finite then AB do exist. It is not difficult to prove 
the following proposition. 

Proposition 3.1. 

(a) A lower or upper triangular matrix A € M.oo(A) is invertible if and only 
if all n-minors of A are invertible. 

(b) An LU -factorization of a matrix A £ M^A) exists if and only if it 
exists for all n-minors of A. Moreover, in this case the LU -factorization 
of A is unique. 

Definition 3.2. Given a matrix A 6 M.oo(A) the shifted matrix s(A) of A 
is the diagonal blocked matrix formed by the identity matrix of size lxl 
and A, that is 

•< A > = C a)- 

To apply an elementary row operation to a matrix A is equivalent to mul- 
tiply A on the left by a lower triangular matrix. For some special matrices 
A e M.oo(A) it is possible to perform the elimination process of Gauss by 
using a periodic sequence of elementary row operations, in the sense that 
there exist a lower triangular matrix To G A4oo(A), that would constitute 
the fundamental part of the periodic sequence, such that the whole periodic 
sequence of elementary row operations is of the form, 

T ; next s(T ); next s 2 (T ); next s 3 (T );--- . 

In other words, this means that the sequence of matrices 

A, T A, s(T )T A, s 2 (T )s(T )T A, s 3 (T )s 2 (T ) S (T )T A, ••• 

converges to an upper triangular matrix. This suggest the following defini- 
tions. 

Definition 3.3. Given a lower triangular matrix To € .Moo (-4) the left 
iterated matrix T L corresponding to To is the infinite (from right to left) 
product 

T L = ---s 3 (To)s 2 (To)s(To)To. 
It is clear that this product converges to a lower triangular matrix. Similarly, 
the right iterated matrix T R corresponding to To is the lower triangular 
matrix given by the infinite (from left to right) product 

To R = T oS (To) S 2 (To)s 3 (To)-- - . 

We say that a matrix A £ MoolA) admits a periodic Gaussian elimination 
process if there exist a lower triangular matrix To € Moo(A) such that TqA 
is upper triangular. In this case To is called a fundamental sequence of the 
Gaussian elimination process of A. 

It is not difficult to prove that Tq and T R are respectively characterized 
by the identities 

T L = s(T L ) T and T R = T s(T R ), 
and it is clear that if To is invertible then T L is also invertible and 
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Now the following proposition is clear. 



Proposition 3.4. An infinite matrix A admits a periodic gaussian elim- 
ination process if and only if A admits an LU- decomposition. In fact, if 
A = LU then the fundamental period is Tq = s(L)L~ 1 and, conversely, if To 
is the fundamental period of A then 

I' 1 = (T ) L and L = (Tq -1 ^ . 

3.2. Relevant infinite matrices. In this subsection we introduce some 
infinite matrices that will be used frequently in what follows and we also 
collect some nice properties that these matrices enjoy. 



The Vandermonde matrix: Vi 



u 



The diagonal matrix formed by the powers of q € A: (D, 



qjij 



The diagonal matrix formed by the factorial numbers: Fij = 5{j (i — 1)!. 
The lower triangular matrix formed by the Pascal numbers: Py = (jZ\)- 
The lower triangular matrix formed by the Stirling numbers of the second 
kind: Sij = S(i,j) for i > j, where S(i,j) = 4j ££ =0 (-l) fc (£) (j - k)\ 



These matrices are, 



V 




D„ 




13 1 
17 6 1 
1 15 25 10 1 



Next we collect some of the properties of these matrices: 

(1) The Vandermonde matrix V has an LZ) [/-factorization. It is given by 
V = PFSK 

(2) P- 1 = D-iPD-i. 

(3) The classical recurrence relations 



o = (;■:» + c- 1 ) 



and 



S(i,j) = S(i-l,j-l)+jS(i-l,j) 



that define, respectively, the Pascal and Stirling numbers correspond to 
the matrix identities P = s(P) To,p and S = s(S) Tq^s where 



To,p 




and 



To,s 




Thus P and S are the left iterated matrices P = Tq P and S 
(4) P and S are also right iterated matrices. Since P = ^ \ \ 1 ; 

/ 1 ■ \ R 
S = Ps(S) we have P=(\\ 1 



1 o,s- 



s(P) and 



and S 
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(5) Combining properties (pQ) and Q we obtain that the Vandermonde ma- 
trix V = PFS 1 admits a periodic Gaussian elimination process with 



fundamental sequence 



i 

l l 
ill 




(6) Let p € M[t] be a polynomial and let pj and Oj be, respectively, the 
coefficients of p in terms of the bases {P} and {{t)j}, that is 

Pit) = Po + Pit + P2t 2 + p 3 t 3 + • • • 

= a (t) + at(t)i + a 2 {t) 2 + a 3 (t) 3 + ... . 

Then the coefficients a* and pj are related by 



, Po 
in . / pi 

P> 



3.3. The matrix representation of the system (jl.3p . Assume that iif 
and E are arbitrary elements of A, y particular we do not assume that 
[H, E] = cE. 

Recall that p = po + p\ t + p 2 t 2 + p% t 3 + . . . satisfies the system 

E n p(H + n)= p(H - n)E n for all n G N, 

if and only if the vector (po>Pi;P2) • • • ) is a solution of the following linear 
system, 

Ep +E Pl (H+l) +E P2 (H+l) 2 ...= Po E + Pl (H-l)E +p 2 (H-l) 2 E .. 
E 2 p +E 2 Pl (H+2)+E 2 p 2 (H+2) 2 . . . = Po E 2 + Pl (H -2)E 2 +p 2 (H -2) 2 E 2 . . 
E 3 p +E 3 Pl (H+3)+E 3 p 2 (H+3) 2 . . . = Po E 3 + Pl {H -3)E 3 +p 2 {H -3) 2 E 3 . . 



or equivalently, that 
f Le LeRh+i LeR h+1 LeR h+1 ... ^ / P0 \ 



L E L E R H+2 L E R H+2 L E R H+2 
L% L 3 e Rh+3 L e R H+3 L E R H+3 



P2 



/ Re ReRh—i ReR^_i ReR?h—i ■■■ \ 
R? E R? E R H -2 R%R%- 2 r2 e r h-2 - 



R% R%R H - 



R E R H-3 



R E R H-3 



( l R H +i R 2 H+1 R 3 H+l 



M = D Ll 



H+l 

1 R H +2 R 2 H +2 R H+2 
1 R H +3 Rr+3 R H+3 



f i Rh-i R'b-i R 



D Re 



■3 

H-l 



1 Rh-2 Rir-t R 3 i 



1 R H . 



H-2 

2 R 3 
H-3 "-H-3 



Now, since R 3 '^ = Y,k>\ {{Z\) (±0 k ~ l R°H [t follows that 
1 1 R H +i R 2 h+i r h+i ■ 1 

1 R H +2 R 2 H +2 R H+2 
1 R H +3 R 2 h+3 R H+3 ' 



^4, 



f 1 Rh-i R-h-i R H-1 
1 Rh-2 R 2 u_i R\ 



1 H?r_, R% 



■po 

VI 

P2 



\: : : : / W V; i i / V i 

where the matrices involved belong to .M 00 (Endk(.M)). Moreover, the above 
identity holds if and only if M I ) =0, where 



H-2 "if-2 

R? 
3 n H-3 



-VD-iPf 
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where 



Hence M 



(S ':)•'•' j 




MqP Rh where 

/ Lr—Rf; 



-Re 



L\-R\ 2L E +2R% 



Le+Re 



V 



T 3 - 

E 



R% 3L 3 R +3R% 



9L%-9R 3 E 27L|+27R| 



/ 



It is also convenient to consider the Pochhammer basis {{t)j} of k[t], 
instead of the basis {t 3 }, to express the system (|1.3j) . If we do this, an 
analysis similar to the previous one shows that p = ao(i)o + +^2(^)2 + 



£ Ai [t] satisfies the system (I1.3D if and only if M' 
M' = (M[-M^)F(P R J, 



0, where 



and the matrices M{, M' 2 and P' x are as follows, 

% j-l)( x )i-j 



M{ = D Ll 



P' 

x 



1111 
12 3 4 
1 3 6 10 
1 4 10 20 



1 



and 



(x) 2 2{x) 1 1 
(x)3 3(a;)2 3(21)1 1 



M' 2 



( 1-1 

1-21 



V 



1-3 3-1 
1-4 6-4 1 



of £j3.2l the matrices M and M' are 



Remark 3.5. According to item 
related by 

M = M'D-^SfD^. 

Thus, the matrix M' could be seen as the result of factoring out from M the 
Stirling numbers that the Vandermonde matrix (appearing in M) contains. 

4. The LU-decomposition of the system (11.31) 



4.1. The L and U factors. We have just seen that M is the matrix corre- 
sponding to system (II. 3D in the canonical basis and M' is the corresponding 
matrix in the Pochhammer basis. We also have 



M 



MoPr h , 



M' = (Mi - M! 2 )F{P' RH ) t : M = M'D-xaiSyD-i. 



Therefore, in order to find the LU-decomposition of the system (|1 .3)) it is 
sufficient to find the LU-decomposition of either M$ or M[ —M 2 . Moreover, 
the entries of Mq and M[ — M 2 are homogeneous polynomials in the variables 
Le and Re- Since an homogeneous polynomial q{xi,x%) in two variables x\ 
and X2 is completely determined by the polynomial q{x) = q(x, 1), we shall 
simplify the notation by substituting the variables 



(4.1) 



Le by x 



and 



Re by 1. 
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Under this transformations, the matrices Mq, M{ and M 2 become 

-1 X— 1 IE+1 

c 2 -l 2i- 2 +2 4x' 2 -4 8z 2 +8 
c 3 -l 3x 3 +3 &r 3 -9 27x 3 +27 



Mi 



o 




and 



1 -2 



V 



1 

1-3 3 
1-4 6 



We point out that the LU-decompositions of the matrices M[ and M 2 , sep- 
arately, are 



M[ = D x PP f and M 2 = P 



( 



\ 



-l l 
i -l 
-l l 
l 



The LU-decomposition of M[ — M 2 is more subtle and it is obtained from 
Theorem 12.21 as follows: 

(1) (M[ — M 2 )F is the matrix corresponding to the system (jl for H = 0, 
when it is expressed in terms of the basis {(t) ,-}. 

(2) In other words, (M[ — M' 2 )F is the matrix corresponding to the system 
Eq° for H = 0, when it is expressed in terms of the basis {{t)j}. 

(3) The definition of the system Eq fc and the Recursion Formula 12.81 imply 
that the system Eq fc is obtained by keeping the first k equations of the 
system Eq fc_1 and by incorporating the equations 



-2R E (e^zl) + R 2 E (eq^ 2 l ) - ad^eq^ 1 ) + R E ad^eq^ 1 ) 



for i > k > 0. Applying the change of variables (14. 1|) these equations 
are transformed as follows, 



eq ?: 



This means that the lower triangular matrix that transforms the system 
Eq fc_1 into the system Eq fc is s fc-1 (T ) with 



(4.2) 



i 

-x-1 



1 

-x-1 1 



(4) Theorem 12.21 implies that the left iterated matrix Tq associated to To 
transforms the system Eq° into the system Eq°°. 

(5) Let 

M[ -M 2 = L'U' 

be the LU-decomposition of M[ — M 2 . According to the previous item, 
we have that V = (T^y 1 = (T^ 1 )^ and U' is of the form U[ - U' 2 
where U[F (respectively, U 2 F) is the matrix corresponding to the left 
(respectively, right) hand side of the system Eq°° for H = (and thus 
c = 0), when the system is expressed in terms of the basis {{t)j}. 
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(6) According to Corollary 12.41 when H = 0, the ith equation of Eq°° is 

i-l 

ec 1 • E [f; 1 ) ad ^ r d i +r ~ 1 p® Er + tr 1 ) ^Y" 1 d^" 1 pm r+i 

r=0 

i-l 
r=0 

Now, replacing ad# by (x — 1) and E by 1 we obtain, 

i-l 
r=0 

i-l 



EC; 1 ) d[ +r ~ 1 P (i-2)( x -iy- r - 1 . 



r=0 



(7) If p = ao(t)o + + ^2(^)2 + • • • then, applying (|2.4|) . we obtain 



e ^ 1; E E (y)(j-i-r + 2) i+r _i(i) J -_ i _ r+ ix(x-l) 



i—r—1 



3 



j>0 r=0 

= E E C; 1 )(i-^-r + 2), +r _ 1 (i-2)^_ r+1 (x-ir r - 1 a r 

j>0 r=0 

(8) Now, the entry of the matrix U[F (respectively, U%F) is the coef- 
ficient multiplying aj—% in the left (respectively, right) hand side of the 
ith equation the system Eq°° given in (7), that is 

(^)« = 77TW E tr') V-i~r+ l) i+ r-i(%-i-r(x - IT 



^ 'r=0 
1 J_i 

(J 



®« = TT^W E V) V-i-r+ l)*tr-i(i - 2) i - 4 - r (x - l)'^- 1 . 



r=0 

We shall now complete the description of the LU-decomposition of the 
system (II. 3D by showing that L', (L') _1 , t/j and [/^ are matrices whose 
entries are given by some particular Jacobi polynomials. The basic facts 
about Jacobi polynomials that we shall need are collected in 

4.2. The entries of L and L are ultraspherical Jacobi polynomials. 

Recall that according to item([5]) of §4.11 we know that L' = (T^ 1 )^ and 
(Z') _1 = Tq. It is easy to see that the inverse of 

(-X-1 1 ■ \ 

X — X— 1 1 
X -x-1 1 ■ I 

is the lower triangular matrix with entries (Tq -1 )^. = x if i > j, that 



is 



/ 1 

/ l+x 



T -l 



1 

1+x+x 2 l+x 1 
l+x+x 2 +x 3 1+x+x 2 l+x 1 
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Theorem 4.1. The entries of Lf and its inverse (L')~ l are given by the 
following ultraspherical Jacobi polynomials. 

(L% = pZ>-\x\ 

((LT 1 ). ={-!)*-' ipitfr)- 

Proof. Let K' and K" be the matrices defined by the right hand sides of the 
above equalities, that is 

K[ 3 = (-1)*-* p^-\x) and = ( l i' ' 3 - p>^(x). 

By the definition of left and right iterated matrices, it is sufficient to prove 
that 

K' = Tq 1 s(K') and K" = s{K") T . 
The first identity holds if and only if 

(-ir j piY(x) = £ xl ~ k x +l _~ l (-i) fe ~ j P^'- k+ \x), 

k=j 

or equivalently, if the following equality holds 

V, % '(*) = g Xl ~'~**\ " 1 (-1)*+^ p- k h - j+l '- k - j+ \x), 

k=0 

which in turn holds if and only if 

pry (f±i) = E (i - (i - ■ v ^ (f±±) . 

fc=0 

Now, this last identity is equivalent to 

p-!f\x) = ((V)' ' / '"' " W ' '■' ■) Pr k - j+1 '- k ' j+ \x), 

k=0 

and this last equality is true, since it is the difference of the two identities 
stated in Lemma 15, II for n = i — j + 1 and a = i. This completes the proof 
of the identity K' = Tq 1 s(K'). 

The second identity holds if and only if 

(-I)*-' \ />;•>'•; = E(-ir* k —\ fct*-Hx) (T ) kj 

k=j 

for i > 2, or equivalently if 

which in turn holds if and only if 
1 ifj (j-i) 

~ i-l r i-j \x-l )^~i-lx-l \x-l J'i-1 (ic-1) 2 r i-j-2 \x-l ) ' 
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or equivalently, if the following equality holds 



A{i-l)j P^Ax) =4*0-1) ^ 



(x)+4ijx fff i _ 1 (x)+iC7+l)(x : '-l) P/l+M+^x). 



Now, this last identity is true since it is (1 + j) times equation (|5.5p with 
n = i — j — 1 and a = j + 1, plus 2(j — 1) times equation (|5.2I) with n = i — j 
and a = j. This completes the proof of the identity if" = s(K") Tq. □ 

Remark 4.2. It is known that the ultraspherical Jacobi polynomials consti- 
tute a 2-parameter family of polynomials that are orthogonal with respect to 
a continuous measure. The above theorem provides the following "discrete 
orthogonality" relationship that involves once many of them, 



(\po 



\ 



1 1,1 2 2,2 

2^1 2^0 

1 1,1 2 2,2 3 3,3 

3^2 3 P 1 3 P 

1 1,1 2 2,2 3 3,3 1 4.4 

^4^3 3P2 4P1 tpo y 



-2,-2 -2,-2 

Pi Po 

-3,-3 -3,-3 -3,-3 
Pi Pi Po 



\P 3 



4,-4 -4,-4 -4,-4 -4,-4 

P 2 Pi Po 



\ 



J 



\ 



4.3. The entries of U{ and \J' 2 are Jacobi polynomials. 

Theorem 4.3. The entries of U[ and U'^ are given by the following Jacobi 
polynomials. 



U' 2 



■ 1 1 <vf' 



Proof. The theorem will be proved if we show that 



IJ 



and 



U' 



I / x+l 



2 I ~ 



■1)3-* (^)- 1 I 



j+2,-2 



(*). 



In fact, according to item 



3-1 



x+l 



'J 



x+l I x^ 
2 \ 2 



2TFT)T E (V) V ~ 1 ~ r + l)i+r-l(0i-i-r 
r=0 

1 )" < ^jrEW(<-'-M0i- 4 -rW) r 



x+l / ag-l 'y-' (-1) J ' 
2 I 2 J (i-i)! 



r=0 



(l-iW-j + r)^ (2fi)' 



r=0 
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u' 



I I x+1 



2 \ — 



n 



Cj-i) 



IT)T E C; 1 ) 0' - i - r + l)i+r-i(i - 2) i _ i _ r 



_1\ r— i+1 



r=0 



= (^r ^ E er) « - rut - (^y 



r=0 
j-i 



= {^r i w E (V) (i - Oro - i + r^- (^)' 

r=0 

= (-ir(¥) ! ^ V2 w- 

This completes the proof of the theorem. 

We have therefore proved the following theorem 
Theorem 4.4. The LU- decomposition of 

x—l x+1 



□ 



X X • \ 

x 2 -l 2a; 2 +2 3s 2 -1 4x 2 

x 3 -l 3:r 3 +3 6x 3 -3 10x 3 +l 

v y 



is M[-M' 2 = L'U' where 

(u% = i-iy-iix-i)*-*- 1 

Corollary 4.5. The LU- decomposition of 

/ x—l x+1 x—l x+1 ■ \ 

x 2 -l 2x 2 4-2 Ax 2 -A 8.-r 2 +S 

M = 



* > i; 

xpjlf(x) -P~lt 2 '~ 2 ( x )) > i< 



J- 



\ 



x 2 -l 2x 2 +2 Ax 2 -A Sx 2 +S ■ 
x z -l 3x 3 +3 &r 3 -9 27:r 3 +27 ■ 



IS 



M = L U 
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/ 1 -1 1-1 1 
1-3 7 -15 
1 -6 25 
1 -10 
1 



where 



(Lo) 
(Uo) 



i > r, 



(-iy-jl( x -i)*-ipjir\ x ), i<j- 



and the explicit numerical matrices are respectively formed by the factorial 
and (up to the minus signs) Stirling numbers of the second kind. 



Proof. We know that 



M = (M[ - MfiFD-^SfD-i 
= L'U'FD^siSfD^. 
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This implies that L' is the L-part of Mq and thus we obtain the expression 
for (Lo)ij. Now we must prove that 

which is easily seen to be equivalent to prove that 




= & 

i i i 

that is 

(x-i) fj pji\ rl (x) + J . (x - 1) pjlt-f 1 ^)^ = x pjif( x )-p~ j 2 2 



3-1 



(x) 



or 



2 P[X 1 (*) + J — 1 P^tY 1 ^)) =(x+l)P~lf(x)-{x-l)PT^-\a 

It is now straightforward to derive this identity by using the formulas (4.5.1) 
and (4.5.4) from [Si]. □ 



5. Appendix: Jacobi and Gegenbauer polynomials 

The Jacobi polynomials are defined for non negative integers n 
arbitrary (rational) numbers a and (5 as follows (see Chapter IV in |Sz] 



P^(x) = ( n \n + a + P + 1) . . . (n + a + (3 + k) 

U ' k=0 ^ ' 

/ i \ k 
fx — 1 

x (a + k + 1) . . . (a + n) 
with the understanding that the general coefficient 

Tl\ 

J (n + a + + 1) . . . (n + a + j3 + k)(a + k + 1) . . . (a + n) 

is equal to (a+l) . . . (a+n) if k = 0, and equal to (n+a+P+1) . . . {2n+a+(3) 
if k = n. The Jacobi polynomials can also be represented as 

P^(x) = (Q + 1) --; (n + a) 2 F a (-n, n + a + /J + 1; a + 1; ^ 
n! \ 2 

where 2-^1 is the hypergeometric function of Gauss. If a = /3 the normalized 
Jacobi polynomials 

r(q + l) r(n + 2a + l) 
r(2a + l) T(n + a + l) n 1 ' 

are called ultraspherical polynomials or Gegenbauer 's polynomials. In what 
follows we shall call ultraspherical Jacobi polynomials to the "unnormalized" 
Gegenbauer's polynomials Pn' a (x). 
Let us consider 

K'^(x) = (x - l) n P^ fe) , 
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it is clear that p" (x) is again a polynomial. These polynomials can be 
expressed in terms of the hypergeometric function of Gauss as follows (see 
(4.22.1) in [Si]), 

Pn (x) = ; 2F1 {—n, —n — a; —2n — a — p;x 

n! 

5.1. Polynomial identities. It is well known that 

(5.1) P«e(x) = {-l) n P^(-x). 

In particular, the ultraspherical Jacobi polynomials are even or odd accord- 
ing as n is even or odd. The ultraspherical Jacobi polynomials also satisfy 
the following identities (see (4.7.14) and (4.7.28) in [Sz]). 

(5.2) (n + 2a)P^ a {x) = 2(n + a)P^ a ~\x) + x(n + a)P^ x {x) 

(5.3) x(n + 2a)P^(x) = 2(n + l)^*" 1 ^) + (n + a)P^{x). 
The difference between these two identities gives, 

(5.4) (l-x)(n + 2a)P^ a (x) 

= -2(n + m:^- 1 (x) + 2(n + a)P^' x (x) - (1 - x)(n + a)P^ (x) . 

On the other hand, multiplying by x equation (|5.2|) and subtracting from it 
equation (|5.3|) yields the following identity, 

(5.5) 2(n+l)P n "~ 1 1 ' a - 1 (x) = 2x(n+a)P^ a -\x) + {x 2 -l){n+a)P:^{x). 
Lemma 5.1. The ultraspherical Jacobi polynomials satisfy, 

k=0 

2j2(^) n - k pr n - k ' a+n - k w = ^ pri*) - p^ix) 

k=0 

Proof. The second identity follows immediately from the first one by using 
()5.ip . then we shall prove the first identity by induction on n. It is clear 
that this identity holds for n = and all a, hence if we assume that it holds 
for n and all a we have, 

n+l 

2 ^ ^l^x^ n+1 - k pa+n+l-k,a+n+l-k ^ 
k=0 

1 _ n 

= l -^2Y J (^) n - k pr n+1 ~ Ka+n+1 -\*) + 2^Ti(*) 

k=0 

1-x n + 2a + 2 ^ a+la+1 , l-x a+lj(X+1 



pr i ' a+i (x) + p:ir +i (x) + 2P^ 1 ( X ) 



2 n+a+l " v/ 2 

n + 1 T ^nTl(^) + W*) + 2P„T 1 (X) 



n + a + 

n + 1 + 2a 
n + 1 + a 
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which completes the proof of the first identity. We point out that in the last 
but one equality we used the identity (|5.4p with a + 1 instead of a. □ 
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